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ABSTRACT 

Asteroid 2008 TC3 was a Near Earth Asteroid that impacted the Earth on 2008 Octo- 
ber 7. Meteorites were produced by the break-up of 2008 TC3 in the high atmosphere 
and at present, about 600 meteorites - caUed Almahata Sitta - coming from 2008 TC3 
have been recovered. A mineralogical study of Almahata Sitta fragments shows that 
the asteroid 2008 TC3 was made of meteorites of different types (ureihtes, H, L, and E 
chondrites). Understanding the origin of this body and how it was put together remain 
a challenge. Here we perform a detailed spectroscopical and dynamical investigation 
to show that the most likely source region of 2008 TC3 is in the inner Main Belt at 
low inclination {i < 8°). We show that asteroids with spectroscopic classes that can 
be associated with the different meteorite types of Almahata Sitta are present in the 
region of the Main Belt that includes the Nysa-Polana family and objects of the Back- 
ground at low inclination. Searching for a possible scenario of formation for 2008 TC3, 
we show that there is little chance that 2008 TC3 was formed by low velocity collisions 
between asteroids of different mineralogies, in the current asteroid belt. It seems more 
likely that the heterogeneous composition of 2008 TC3 was a inherited from a time 
when the asteroid belt was in a different dynamical state, most likely in the very early 
Solar System. Because ureilites are fragments of a large, thermally metamorphosed as- 
teroid, this suggests that the phases of collisional erosion (the break-up of the ureilite 
parent-body) and collisional accretion (the formation of the parent body of 2008 TC3) 
overlapped for some time in the primordial asteroid belt. 
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1 INTRODUCTION 

Meteorites are a partial sample of asteroids that survive the 
passage through the Earth's atmosphere. The identification 
of the source regions of the different type of meteorites is 
essential to be able to link the mineralogical properties of 
meteorites to the parent asteroids and, consequently, to ad- 
dress the mineralogical evolution in the asteroid belt. How- 
ever, this is not an easy task and only some of these links 
could be established: for instance, the group of howardites, 
eucrites, and diogenites (HEDs) meteorites ar e thought to 
come from the Vesta family of asteroids (e.g. iBinzel fc Xul 
I1993I ): more speculativel y, L ordinary chondrit es could come 
from the Gefion family l|Nesvornv et al.ll2009l 'l. while aster- 
oids of the the Fl ora family bear spect ral similarities with 
the LL chondrites jVernazza et ahlboOSl ). However, the par- 
ent bodies of most meteorite types, if still intact, are un- 
known. 

The discovery and spectroscopic observation of the 
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near-Earth asteroid (NEA) 2008 TC3 (henceforth TC3) 
20 hours before it impacted the Earth's high atmosphere, 
and the subsequent recovery of meteorites (called Almahata 
Sitta) - clearly com ing from this body - was a major re- 
sult in this respect (jJenniskens et all |2009| : IShaddad et al.l 
It allowed a direct link between an asteroid and 
meteorites to be established for the first time: the aster- 
oid was classified as belongin g to the spectroscopic F-class 
(in the Tholen classifi cation, Tholenl Il 984l or B-class (in 
the Bus classification, 'Bus fc Binzcl 200!^) on the basis of 
the fiat shape of its reflectance spectrum in the region be- 
tween 500 and 1000 nm. Moreover, among the 47 meteorites 
initially recovered, it was observed that the visible spec- 
trum of the fragment #7 matches the telescopic spectrum 
of TC3 obtained before the im pact with the Earth's atmo- 
sphere jjenniskens et al] 20091'). Fragment #7 is an achon- 
drite polymict ureilite ( Jenniskens et al"]|2009f ). Ureilites are 
in a group of achondritic meteo rits that exhibit bot h primi- 
tive and evolved characteristics (jCloutis et al.l|201o| ): in par- 
ticular, they are characterized by olivine and pyroxene-rich 
clasts among carboneous material (mainly graphite); flne- 



2 J. Gayon-Markt et al. 



grained graphite is also present, which lower the albedo 
of the meteorites (about 10-12%; iHiroi et al.1 12OI0I '). Ure- 
ilites were initially thought to derive from S-class asteroids 
(see for instance Gaffcy ct al. 1993i). Howeve r, because of its 
spectral properties, Ijenniskens et all l|2009l ) propose a link 
with B-class asteroids according to the Bus classification or 
F-class in the Tholen classification. This is more plausible 
than a link with S-class asteroids, given the low albedos of 
ureilite meteorites, consistent with B or _F-class asteroids 
but not with the S'-class. 

It is worth to note that the F-class can be distinguished 
from the B-class from a much weaker UV drop-off in the 
spectra of the former compared to the latter and also be- 
cause B-class asteroids show a moderately higher average 
albedo than F-class bodies. However, in the Bus classifi- 
cation these two classes are merged in a unique class (the 
B-class). This is because the IBus fc B inzcl (200^) spectral 
classification is based on spectra acquired with CCD spec- 
trographs, which - in general - do not observe far enough 
in the UV to observe the above-mentioned drop-off feature. 
We will refer in this work to B-class asteroids including both 
Tholen B- and F-classes, and Bus B-class. 

Interestingly, mineralogical studies of Almahata Sitta 
show that, among the ~ 600 fragments recovered 
l|Shaddad et al.ll2010l 'l. about 70 - 80% are ureilites, while 
the remaining 20 — 30% are enstatite chondrites, H and L 
ordinary chondrites. More specifically. iBischoff et all (|2010l ) 
show that, among a subsample of 40 deeply studied mete- 
orites from Almahata Sitta, 23 fragments are achondritic 
ureilites and 17 have chondritic litologies with 14 of them 
corresponding to enstatite chondrites, 2 to H ordinary chon- 
drite s and 1 to a new type of chondrite (see lHorstmann et al.l 
I2OIOI for more details). Although small clasts of differ- 
ent types arc quite common in brecciated meteorites (see 
iMeib om & Clark 1999, for a review) and carboneous mate- 
rial is found in some HED meteorites, it is the first time that 
meteorites of very different mineralogies (i.e primitive frag- 
ments with achondrite polymict ureilites and evolved ones 
such as ordinary chondrites or enstatite chondrites) are as- 
sociated to the same fall. This led to make the hypothesis 
that TC3 was an asteroid made of blocks of different min- 
eralogies, held together very loosely (give n the explosion of 
the body at the anomalously high altitude ; Jenniskens et al.l 
I2OO9I : iBischoff et allbOlOl : IShaddad et al.ll2010l ). 

Tracing back TC3 to its source region in the asteroid 
Main Belt would allow us to understand the origin of the 
Almahata Sitta meteorites and how TC3 was put together 
by loosely assembling meteorites of different mineralogies. 
Establishing this link would also be fundamental to shed 
light on the source region of ureilites, that albeit rare, are 
the forth major class of primitive meteorites recovered on 
Earth after the CV, CI and CO carbonaceous chondrites. 

In their att empt to find the source region of TC3 and 
Almahata Sitta, Ijenniskens et al] (|2010t ) selected all B-class 
asteroids, according to the Bus classification, and objects 
of the Tholen F and B classes and searched for spectra 
similarities with TC3 and Almahata Sitta. As a result of 
their study, they showed spectral similarities between TC3 
and ungrouped asteroids as well as several dynamical as- 
teroid groups (or families) as possible origins for the TC3 
asteroid, namely Polana (2.4 AU, 3°), Hoffmeister (2.8 AU, 
4.5°), Pallas (2.8 AU, 33°), Themis (3.15 AU, 1.5°), and 



Theobalda (3.2 AU, 14 °). Later, from dynamical grounds, 
Ijenniskens et all (|2010l ) concluded that asteroids from the 
inner asteroid belt (i.e. with orbital semimajor axis a < 
2.5AU) are the likely parent bodies of TC3. This reduces 
the choice to dispersed B-class asteroids in the inner main 
belt and the Polana asteroid group. In Section (21 we revisit 
this issue studying the possible dynamical source regions for 
TC3. 

We recall here that the Polana group is part of a clus- 
ter of asteroids known as the Nysa-Polana family l|Nesvornvl 
l2010t ). which is located in the inner main belt, between the 
!/6 secular resonance (at heliocentric distance « 2.1 AU) and 
the 3:1 mean motion resonance with Jupiter (at heliocentric 
distance of 2.5 AU). This family has a co mplex - twofold - 
structure in orbital proper element space (|Nesvornvll2010l ). 
suggesting that it is the outcome of at least two independent 
break-up events in the same orbit al region. From the few 
spectral data available at the time, ICellino et al.l (|200ll ) ar- 
gued that the Nysa-Polana family contains asteroids of three 
spectral classes. The first class is that of B- class objects, 
like asteroid (142) Polana itself - note that [Cellino et ah] 
(2001) uses the F-class classification from thejTholcn (lOSj) 
taxonomy; the second class is the S-class, with the largest 
member being identified as the asteroid (878) Mildred; the 
third class is that of X-class objects, such as the asteroid 
(44) Nysa. In this manuscript, we revisit this result using 
a much wider dataset of spectro-photometric data provided 
by the Moving Objects C atalog of the Sloan Digital Survey 
(SDSS, E vezic et al.ll2002l ). w hich is analyzed he re using a 
new classification algorithm (|Michel et al.l [2OO5I . described 
in Section |3J developed for the Gaia space mission of the 
European Space Agency. 

A detailed study of the mineralogy of the Nysa-Polana 
family is of great importance also for bett er understand- 
ing th e origin of other NEOs. In partcular, ICampins et al.l 
(|2010h claimed that the asteroid (101955) 1999 RQ36, tar- 
get of the sample return mission OSIRIS-REx (approved 
in the program New Frontiers of NASA), was delivered 
to near-Earth space via the secular resonance from the 
Polana group. Moreover, the binary asteroid (175706) 1996 
FG3, primary target of the sample return mission Marco 
Polo-R, under study at the European Space Mission (ESA), 
might have formed within the Polana group and delivered 
to the near-Earth space via the overlapping Jupiter 7:2 and 
Mars 5:9 mean m otion resonances rather than the (see 
IWalsh et al1l2012l ). 

As a consequence, in Sections U and [3 we perform a 
spectroscopic analysis using the SDSS data of the asteroids 
of the Nysa-Polana family as well as dispersed asteroids of 
the inner Main Belt (called objects of the background) in 
order to find spectral matches with TC3 and Almahata Sitta. 

Finally, in Section [6l we investigate a possible forma- 
tion scenario for the TC3 asteroid as a rubble pile of rocks 
of different mineralogical types, which is based on the pecu- 
liar low inclination of the Nysa-Polana family and dispersed 
asteroids. 
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2 DYNAMICAL HISTORY AND MAIN-BELT 
ORIGIN OF TC3 

The large majority of NEAs are fragments generated by the 
coUisional disruption of larger asteroids of the main belt; 
said fragments drift in orbital semimajor axis by the so- 
called Yarkovsky effect until they reach regions of orbital in- 
stabilities (mean motion resonances with Jupiter and secular 
resonances) which, by enhancing their orbi tal eccentricities, 
deliver them to the near-Earth space (see iMorbidelli et alJ 
[2002 . for a review). 

According to model by iBottke et al.l (|2002l ) , there are 
five main intermediate sources of NEAs: 

1) the Ufj secular resonance, which marks the inner edge of 
the asteroid belt and occurs when the precession frequency 
of the longitude of perihelion of an asteroid is equal to that 
of Saturn; 

2) the 3:1 mean motion resonance located at a ~ 2.5AU, 
where the orbital period of an asteroid is 1/3 of that of 
Jupiter; 

3) Mars-Crossing asteroids, defined as objects which are 
not NEAs (i.e. their perihelion distance q is larger than 1.3 
AU), but whose semimajor axis evolves in a random- walk 
fashion as a result of close and distant encounters with Mars. 

4) the outer belt population, whose eccentricities can 
increase up to planet-crossing values due to a net- 
work of high order orbital resonances with Jupiter 
and three-body resonance of type asteroid- Jupiter-Saturn 
ijMorbideUi fc Nesv ornv"l999h 

5) dormant Jupiter family comets. 

The orbital elements of TC3 before impact are not 
known very precisely. Here are different estimates: a — 
1.29 AU, e = 0.299, and i = 2.441° (NEODyS website: 
^http://newton. dm. unipi.it/neo dys); a = 1.308 AU, e = 
0.312, and i = 2.542° ( Jcnniskc ns et al.ll2010l '). For each of 
these orbits we computed, usin g the original orbital evolu- 
tion files of iBottke et all l|2002l ). which are the most likely 
intermediate sources among those enumerated above. We 
found that TC3 has a probability of 63-66% to come from the 
ve source and 34-37% to come from the Mars crosser popu- 
lation. A pparently, none of th e simulations for the 3/1 res- 
onance in IBottke et al.l l|2002l ) produce objects within ±.05 
AU in a, ±.05 in e and ±2. 5° in i of TC3. This co ntrasts 
with the claim reported in Ijenniskens et all (|2010l ). that 
TC3 has a 20% probabilit y to come from the 3/1 resonance 
in the lBottke et al.l (|2002l ) model. We s uspect that th e re ha s 
been an error in the manipulation of the lBottke et al.l (|2002l ) 
model in that work. 

Of the Mars cro ssers that can produce the TC3 orbit in 
iBottke et all (|2002l ) simulations, none has semimajor axis 
larger than 2.5 AU. From this, and remembering that the 
resonance lays at the inner edge of the b elt, we conclude - 
in agreement with Ijenniskens et al.l (|201Cll ) - that TC3 most 
likely comes fro m the inner asteroid bel t, inside of 2.5 AU. 

As noted in ljenniskens et "al] (|201Cll ). there are two pop- 
ulations of asteroids in the inner belt with spectra broadly 
consistent with that of TC3: the Polana members, with an 
orbital inclination of 2-3 degrees, and a population of dis- 
persed B-class asteroids, with inclinations ranging up to 15 
degrees. These objects are too dispersed to belong to a rel- 
atively young coUisional family, but their broad cluster in 



orbital space suggests that they might belong to an old col- 
lisional family, dynamically dispersed probably during the 
phase when the orbits of the giant planets changed substan- 
tially, about 4 Gyr ago. 

To determine whether it is more likely that TC3 comes 
from the Nysa-Polana family or from the population of dis- 
persed B-class asteroids, we turn, once again, to the original 
simulations in Bottkc ct al. (2002). In that work, the initial 
distribution of asteroids in the i/a was uniform in inclina- 
tion. In particular, 50% of the initial conditions had initial 
inclination larger than 8 degrees. However, 93% of the parti- 
cles which reproduced the orbit of TC3 at some time during 
their evolution, have initial inclination z < 8°. For the Mars 
crossers, all those reproducing the orbit of TC3 have initial 
i < 8° . This suggests that the most likely source of TC3 is 
at low inclination, consistent with the Nysa-Polana family 
and with the dispersed B-class asteroids with i < 8°. This 
is not surprising, given that the inclination of TC3 is 2.3- 
2.5 degrees. However, the dynamics rules out a TC3 origin 
from the dispersed B-class asteroids at higher inclination 
(i.e. i> 8°). 

Finally, ICampins et al.l (|2010l ) have shown that the 
Polana group can easily deliver small enough fragments 
in the va resonance. In fact, the extrapolated Yarkovsky- 
induced semimajor axis distribution of Polana group mem- 
bers predicts that asteroids fainter than H ~ 18.5 can reach 
the border of the i/q resonance, which is at 2.15 AU for a 
Polana-like inclination. Assuming a Polana-like geometric 
visible albedo of 0.05, H = 18.5 translates into a diameter 
D ~ 2 km, which is much larger then the size of TC3 . 

For all these reasons we conclude that, from the dy- 
namical viewpoint, the Nysa-Polana family and dispersed 
B-class asteroids with i < 8° are the most likely sources of 
TC3. Thus, we will focus on the Nysa-Polana family and the 
Background at low inclination in the next sections. 



3 SPECTROSCOPIC ANALYSIS METHOD 

3.1 Asteroid spectral classification algorithm 

In order to search for spectral groups in the Nysa-Polana as- 
teroid family or in the inner main belt background at low in- 
clination, we used the unsupervised classification algorithm 
that will be adopted for the classifica tion of asteroid sp ec- 
tra from the multiband photometers (|Jordi et al.ll2O10l ) on 
board the Gaia space mission of the European Space Agency 
(ESA). 

While Gaia is mainly devoted to the observations of 10^ 
stars, it is expected that this mission will also observe more 
than 250, 000 asteroids multip le times over 5 years (mission 
lifetime) (iMignard et al.ir2007l ). Spectral classification of as- 
teroids will be performed by an u nsupervised c l assific ation 
algorithm based on t he works of I Michel et al. 1 l|2005l ) and 
iGalluccio et al.l l|2008l ). An unsupervised approach has the 
advantage that no a priori information is taken into account 
to build the spectral groups. 

The algorithm is based on a method for partitioning a 
set V of A'^ data points {V G R^) where L in our case is 
the number of spectral bands into K non-overlapping clus- 
ters (or groups of data points) such that (a) the inter-cluster 
variance is maximized and (b) the intra-cluster variance is 
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Figure 1. Minimum distance between spectra computed by the 
MST method as a function of the spectra index. The cut off used 
to obtain the B, X, and S groups was determined such as already 
classified asteroids belong to their appropriate group. 



minimized. Each spectrum is considered as a vertex of tree. 
A tree is a graph that is connected (i.e. every vertex is 
connected to at least another vertex) and acyclic (i.e. the 
graph has no loops). The spanning tree (MST, i.e. the tree 
passing through each vertex of the set) with the minimal 
length is calculated u sing the so-called Prim's algorithm (see 
iGalluccio et al.ll2009l . and references therein). The length of 
each edge of the tree (i.e. the distance between two spec- 
tra) is determined using the KuUback-Leibler metric: let 
Vi = {vi\, . . . , ViL} the feature vector corresponding to a re- 
flectance spectrum (L being the number of wavelengths); at 
a given wavelength A-,, each spectrum is associated with a 



J Vij , which 



(positive) normalized quantity: vij = Vij/^ 
can be interpreted as the probability distribution that a 
certain amount of information has been measured around 
the wavelength A. To measure the similarity between two 
probability density functions, we compute the symmetrized 
Kullback-Leibler divergence: 



dKL{Vi\\vk) = y^iVij ~Vkj)\og^ 



Vkj 



The identification of the cluster is performed by first 
computing, at each step of the MST construction, the length 
of the newly connected edge and then by identifying valleys 
in the curve obtained by plotting the MST edge length as 
a function of the iteration of the construction (see Figure 
The valleys in this curve identify the number and the 
posi tion of high density region of points, i.e. the clusters 
(see IGalluccio et al.l |2009| . for a more thorough description 
of the algorithm). 



3.2 Spectral data 

In order to increase the sample statistics compared to 
the relatively sma ll number of asteroids observed by 
ICellino et aP ()200lh and identify meaningful spectral groups 
within the Nysa-Polana family, we used our spectral clas- 
sification algorithm on the observations contained in the 



SDSS MOC4 (jlvezic et al.ll2002l .lhttp://www.sdss.org). The 

MOC4 contains magnitudes of 104449 Main Belt asteroids. 
Each asteroid was observed, in general, at multiple epochs 
over 5 spectral bands in the visible light, namely u' , g' , r' , 
i' , z at the following central wavelengths (A) of 354, 477, 
623, 763, 913 nm, respectively. 

To obtain the surface reflectance of the asteroids, we 
needed to remove the solar contribution (uq, Qq, tq, i©, 
zq) to the observed magnitudes. The solar contribution 
is calculated from transformation equations between the 
SDSS u'g'r'i'z' magnitudes and the usual UBVRcIc sys- 
tem. We find: u© = 6.55, = 5.12, tq = 4.68, — 
4.57, Zq = 4.54. (See the SDSS website for more details: 
http://www.sdss.org). In order to compare asteroid spec- 
tral reflectances and be able to classify them, we choosed to 
normalize each spectrum relatively to the r' band (A^' = 623 
nm). Finally, we obtained the following asteroid reflectance 
for the A„/ wavelength: 



= 10 
with C,,i 



-0.4(C„ 



■ u — uq and 



the reflectance colors at A^/ and A^/ respectively. The same 
computations were performed for the Xgi , X^i , and A^/ wave- 
lengths. 

We searched in the M0C4 database for all asteroids be- 
longing to the Nysa-Polana family and for dispersed objects. 
Concerning the family members idcntiflcatio n, we used the 
definition of dynamical families of Nesvorny l (j201Gl ). Among 
the Nysa-Polana family, 4134 objects have been observed at 
least once by SDSS (as a comparison, the number of spec- 
tra obtained in the visible light or near-infra red, for the 
Nysa-Polana family, in the SMASS, EGAS and 52-Golor cat- 
alogs are : 15, 13, and 2, respectively). For those asteroids 
with more than one observation, we calculated the weighted 
mean refiectance by averaging the reflectancies derived at 
each epoch weighing each point with its respective signal- 
to-noise ratio. 

However, likely due to non-photometric conditions some 
of the observed objects have high uncertainties in the mea- 
sured magnitudes. In the visible light, the separation of 
asteroid spectroscopic group is based on the overall slope 
of their reflectance and in the presence (or absence) and 
strength of the 1 \im absorption feature. In order to avoid the 
superposition of spectral groups we selected asteroid with 
the less noisy observations. Namely, we rejected observations 
with a relative uncertainty > 10% on the in-band photomet- 
ric flux derived from the MOG4 magnitudes. In the end, the 
taxonomic classiflcation of the Nysa-Polana dynamical fam- 
ily was performed over 579 objects and over 2828 objects for 
the Background at low inclination. 



4 SPECTRAL ANALYSIS OF THE 
NYSA-POLANA FAMILY 

4.1 Application of the spectral classification 
method 

Figure [1] shows the MST length as a function of the num- 
ber of iterations (i.e. addition of spectra) in the case of the 
SDSS-MOG4 spectrophotometric data of the 579 asteroids 
of the Nysa-Polana family. The flgure clearly shows three 
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Figure 2. Panel (a), (b), and (c) represent the spectra for each 
group with their mean spectrum. A comparison with several spec- 
tra classes coming from the Bus taxonomy is also performed: 
Group 1 with S-class spectrum, Group 2 with B-class spectrum 
and Group 3 with X-subclasses. We also find a good agreement 
with several Almahata Sitta fragments. Spectral differences with 
fragment #25 (H-chondrite) appearing in panel (a) are likely due 
to space weathering. Fragments #4 and #7 in panel (b) are ure- 
ilites. 



main valleys corresponding to three groups composed of 118, 
13, and 378 asteroids, respectively. The remaining objects 
(70) are considered as unclassified asteroids. This is due to 
their large value of the MST length which means no spectral 
similarity among the family is found for these 70 objects. 

Figure [2] shows asteroid spectra of each group and 
their average spectrum. Error bars are the 1-a stan- 



Figure 3. Distribution of S-class (■), B-class •), and X-class 
(a) objects of the Nysa-Polana family as a function of a* and 
i' — z' (calculated from SDSS magnitudes). Ellipses represent the 
location of S , B, an d X regions. The color palette is the same of 
iParker et al] ||2008| '). Although Polana (B), Nysa (X), and Herta 
(X) have not been observed by SDSS, we have calculated their 
a* and i' — z' values from the EGAS catalog and located them in 
this plot. 



dard deviation. We also compared the spectra of each 
spectral group of asteroids with the mean spectrum of 
different taxonomic classes from the Bus classification 
(http://smass.mit.edu/busdemeoclass.html). In each panel, 
we only plot the Bus mean spectra which present the closer 
similarity with the mean spectra of each group: 

(i) For the first group (top panel of Figure[2l 378 objects) , 
we find that all spectra resemble the 5-class mean spectra 
of Bus. 

(ii) In the middle panel of Figure[2] the 118 spectra of the 
second group fits with B-class asteroids (such as 142 Polana, 
for instance). 

(iii) The bottom panel shows the third more populous 
group containing 13 asteroids. Comparing the spectra of 
this group and the mean spectrum of the X-class, Xc-class, 
Xfe-class, and Xe-class, we find that the third group of 13 
objects corresponds to asteroids of the X-complex and its 
subclasses. 

We note that groups (1) and (2) are consistent wit h 
the two main spectral classes found bv lCellino et al.l l|200l[ ). 
Moreover, the same authors had already identified the pres- 
ence of asteroids with X-class spectral reflectances. They 
also noted that two of the largest asteroids, (44) Nysa and 
(135) Herta have spectra consistent with an X-class. As a 
consequence, we confirm their results using the SDSS MOC4 
large sample of spectra. 

In Figure [3] we plot the asteroids of the three spectro- 
scopic groups of the Nysa-Polana family as a function of a* 
and i' - z' (with o* = 0.89(5' -r') + 0.45(r' - i') - 0.57) . We 
show that the S-class group is well defined with a* G [0; 0.25] 
and i' ~ z' £ [—0.15; 0.08]. The B-class region stretches from 

0.2 to ~ -0.06 for a* and from -0.1 to 0.14 for i' - z 

and we find for the X-class region, a* £ [—0.1;— 0.02] and 
i' - z g [-0.02; 0.0 8 ]. Th e color palette used in Fig. O is 
that of iParker et al.l (|2008l ). Blue dots correspond to aster- 
oids with neutral or slightly blue spectra (B-class) , asteroids 
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with neutral to slightly red spectra (X-class) are displayed 
in the plot with purple colors, while S-class asteroids are red 
to yellow. 



4.2 Comparison with Almahata Sitta meteorites 

It is very interesting to note that the three asteroid spectro- 
scopic groups that we found in the Nysa-Polana family are 
likely to be the analogs of the different meteorite mineralo- 
gies found in Almahata Sitta. To strengthen this argument, 
we also compare - in Figure [5]- the mean spectrum of each 
asteroid spectroscopic class of the Nysa-Polan a family (S, 
B, X ; red curves) with published spectra (from lHiroi et all 
|2010D of some fragments of Almahata Sitta. We note that 
at the time of writing, spectra of E-chondrites from Alma- 
hata Sitta are not yet publicly available. So, comparisons 
between the classes of asteroids in the Nysa-Polana family 
and fragments of Almahata Sitta were performed for the S 
and B classes only. However, comparison of sp ectra of other 
E-chondrites (from the RELAB database from lGaffevlll976l ') 
with the average spectrum of the the Nysa-Polana X-class 
shows a good agreement. 

Concerning the B class, we find that the mean spectrum 
of the Nysa-Polana B-class asteroids matches the spectra of 
the fragments #4 and #7 (ureilites) of Almahata Sitta. It 
is important to note that this spectral match was obtained 
only considering the visible wavelengths between about 350 
and 900 nm. We remind that, from spectro scopic and albedo 
similarities, Ijenniskens et all l|2009l . |2010| ) have proposed a 
link between the B-class (or the F-class) and ureilite me- 
teorites. An important caveat is that the link between B- 
class asteroids and the Almah ata Sitta ureilites proposed by 
Ijenniskens et al" I l|2009l . I2OIOI ) is based on the noisy spec- 
trum of TC3. Moreover, B-class asteroids are more com- 
monly assoc iated with carbon aceous chondrites by several 
studies (see I Clark et al.l I2OI0I . for a review about B-class 
objects, their meteorite analogues and their composition). 
Ho wever, a recent spect roscopic survey of B-class asteroids 
bv lDe Leon et al.l ((2012') shows that the ensemble of the re- 
flectance spectra of the 45 B-class asteroids analyzed in their 
work have a continuous shape variation in the range between 
500 and 2500 nm, from a monotonic negativ e (blue) slope 
to a positive (red) slope. iDe Leon et al.l l|2012l ) apply a clus- 
tering technique to reduce the ensemble of the spectra to 
6 optimized averaged spectra or "groups" . Interestingly the 
RELAB spectrum of the fragment #4 of Almahata Sitta 
shows a good match with the group #3 of iDe Leon et al.l 
l|2012l ) in the region between 500 and 2500 nm (De Leon, J. 
private communication September, 2011). 

Concerning the S class, we show, in Fig. [2] that frag- 
ment #25 (an ordinary chondrite) is rather close to the 
S-class mean spectrum, as expected, but not in very good 
agreement (see e.g. the review of Chapman 2004, for spectro- 
scopical links between ordinary chondrites and S-class aster- 
oids). One possible explanation for this spectral mismatch 
is the space weathering of S-class asteroids. We discuss this 
possibility in the following subsection. 
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Figure 4. Spectral slope as a function of the exposure to space 
weathering of S-class objects of the Nysa-Polana family (black 
dots and black line) and S-class objects of the background at 
low inclination (circles and dashed line). As a comparison, the 
spectral slope of H-chondrites which are thought to come from 
S-class asteroids is of about 0.1 — 0.2. 



4.3 Space weathering of S-class members of the 
Nysa-Polana family 

Space weathering is a physical process caused by cosmic 
rays, collisions, ion bombardment, that alters physical and 
spectral properties of the surface of atmosphere-less plane- 
tary bodies. More particularly, due to space weathering, S- 
class asteroids become darker (their albedo is reduced), red- 
der (the reflectance increases with increasing wavelength), 
and the depth of absorption bands are reduced. As a con- 
sequence, the spectral slope chan ges due to this process 
dClark et al.ll2002l: [Chapmanll2004 ). 

In Marchi et al.l ( 20061 ). a strong relation was found be- 
tween the spectral slope of S-class asteroids and their "ex- 
posure" to space weathering. The exposure of an asteroid 
to space weathering corresponds to time integral the flux of 
solar ions that the body receives along its orbit. The expo- 
sure to space weathering depends on the age of the aster- 
oid and its average distance from the Sun. Therefore, in or- 
der to know if space weathering affected S-class asteroids of 
the Nysa-Polana family, we calculate t he exposure to space 
wea thering of these bodies as shown in iMarchi et al.l (I2OO6I ) 
and lPaohcchi et al.1 (|2007l ). 

We compute the slope a of a SDSS spectrophotometric 
data by fitting a linear equation normalized to 1 at \r (623 
nm) to the datapoints. To calculate the exposure, we need to 
estimate the average coUisional age of Main Belt asteroids 
(Tmba), which can be constrained by the age of the Belt 
{tLHB ~ 4 Gyr) and by coUisional times (tcou)'- 
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According to lBottke et al.l l|2005h . the coUisional lifetime can 
be estimated as a function of asteroid diameter (using SDSS 
H absolute magnitudes and assuming an albedo of py — 0.2 
for all S-class objects). In the end, we have: 
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Figure 5. Plot 1: Cumulative number of asteroid A'^(< H) as a function of tlie H absolute magnitude for B-class asteroids of the 
Nysa-Polana family (black curve) and the Background at low inclination (i < 8°; dashed curve). Plot 2-4: Distribution in eccentricity 
(plot 2), inclination (plot 3) and H absolute magnitude (plot 4) as a function of semimajor axis. Black points correspond to Polana 
B-cIass asteroids and black circles to B-class dispersed objects from the Background at low inclination. Plot 4 shows the //-dependent 
semimajor axis distribution induced by the Yarkovsky effect that best fits the boundaries of the observed distribution for the Polana 
group (black curves) and for the Background at low inclination (dashed curve). 



with a and e the asteroid proper semimajor axis and eccen- 
tricity, respectively. 

In Figure |4l we plot both the spectral slope of Nysa- 
Polana S-class objects as a function of the exposure to space 
weathering and the best linear fit to this distribution. The 
two-tailed probability for the linear correlation is lower than 
0.001% which means the linear correlation is significant 
This shows that the spectral slope increases with the expo- 
sure. 

Semimajor axes and eccentricities do not vary much 
within the Nysa-Polana family. As a consequence, the expo- 
sure variation mainly comes from asteroid diameter, through 
the computation of Tmba- Hence, plotting the spectral slope 
as a function of the exposure, or Tmba, or asteroid diam- 
eters (d) provide quite similar figures (not shown here) - 
where the exposure increases with time or diameter - but 
with the following abscissa ranges: 0.5 < Tmba < 2.5 Gyr 
or 1 < d < 6 km. 

The spectral slope increase with exposure or asteroid 
age is then the proof that space weathering occurs for S- 
class asteroids of the Nysa-Polana family. We also note that 
the spectral slope of the fragment #25 of Almahata Sitta 
has the value of 0.03. 

As a consequence, the spectral mismatch between as- 
teroids of the S-component of the Nysa-Polana family and 
the ordinary chondrites of Almahata Sitta - as observed in 
Figure [2]- can be explained by space weathering. 



5 ANALYSIS OF ASTEROIDS OF THE 

BACKGROUND AT LOW INCLINATION 

As mentioned in Sections [1] and [2] the other possible ori- 
gin source for TC3, from a dynamical point of view, is the 
population of dispersed asteroids at low inclination. In the 
present section, we perform the same spectroscopic study 

^ A linear fit is considered significant whenever the two-tailes 
probability is lower then 5%. 



as for the Nysa-Polana family (Section |4]) but using ob- 
jects of the Background located in the Inner Main Belt 
(2.1 < a < 2.5 AU) with a proper inclination lower than 8°. 
Selecting the best observations of the SDSS MOC4 catalog 
(i.e. observations with a relative uncertainty > 10% on the 
in-band photometric fiux, corresponding to 2828 objects), 
we find 4 large groups of asteroids corresponding to the fol- 
lowing spectral classes: B, C, X, S, as well as small clusters 
of Q, V, and A-class asteroids. In addition, for the B, X, 
and S groups of the Background, we obtain mean spectra 
very similar to those of the Nysa-Polana family (not shown 
here). We also note that space weathering is also found for 
S-class asteroids of the Background (see Fig. |4]). As a con- 
sequence, the three ingredients required to form TC3 (S- 
class asteroid/H-chondrite; B-class asteroid/ureilite; X-class 
asteroid/enstatite chondrite) are both found in the Nysa- 
Polana family and in the Background at low inclination. 

In our attempt to definitively conclude on the origin re- 
gion of TC3, we have compared the distribution of B-class 
asteroids for the two source regions. In order to get a homo- 
geneous selection of B-class asteroids both from the Nysa- 
Polana family and the Background at low inclination, we 
have selected objects such as their spectra are similar to the 
B-class mean spectrum of the Nysa-Polana family, within 
the 1-a error bar. Although this selection reduces the size 
of the original B-class groups, we limit the overlap of classes 
(B and X classes can overlap; see Fig. [2]). 

From this selection, we plot in Fig. [S] the cumulative 
number of asteroids as a function of H absolute magnitude 
as well as the distribution in eccentricity, inclination and H 
magnitude as a function of semimajor axis. We find more 
B-class asteroids from the Background at low inclination 
(dashed curve in the first plot of Fig. [5]) than from the Polana 
group (black curve). As a consequence, the Background at 
low inclination seems to dominate the Polana group, even if, 
around magnitude 17, the slope of the Polana group curve 
looks a little steeper than the curve for the Background. 
In the other plots in orbital elements, we can see that the 
two regions overlap in some points and that B-class objects 
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of the Background are much more dispersed compared to 
Polana aste roids. 

As in ICampins et"all \2Q1(\ ). we also represent the 
absolute magnitude {H) a a function of semimajor axis, 
bo th for the Polana gro up and the Background (the plot 
of ICampins^ et al.l l2010l . was done for the whole Nysa- 
Polana family). The plot appears to be V-shaped which 
is a feature known to be associated both with the size- 
dependent ejection velocity field and with the drift in the 
proper semimajor axis in duced by the Yarkovsky effect (e.g. 
IVokrouhhckv et ahllioo^ . In the plot, the black curve shows 
the H-dependent semimajor axis distribution induced by the 
Yarkovsky effect that best fits the boundaries of the observed 
distribution for the Polana B-class asteroids. Objects below 
this curve are expected to be interlopers and may belong to 
the Background. The extrapolated Yarkovsky-induced dis- 
tribution predicts that the Polana group should reach the 
outer edge of the resonance for objects with H ~ 18, 
which for a Polana-like albedo of p„ = 5% translates into a 
diameter D ~ 2 km. This means that objects smaller than 
2 km, such as TC3, can easily escape the Polana group and 
the Inner Main Belt. 

Concerning the Background at low inclination, differ- 
ent size limits can be computed due to a large variation of 
proper orbital inclinations. For an inclination of ~ 0°, the 
resonance boundary is found at a ~ 2.1 AU, which gives 
us an H magnitude of ~ 16 and an asteroid diameter limit 
of d ~ 4 km. For the highest inclination (i ~ 8°), we obtain 
H ~ 14.5 at a ~ 2.2 AU, which is equivalent to a B-class 
asteroid (p^ = 5%) of 8 km. As a result, the Background at 
low inclination could have also delivered TC3 through the 
v% resonance. 

Because the Background at low inclination is not a dy- 
namical family, its V-shaped structure in H{a) was not espe- 
cially expected. We then think that the Background at low 
inclination could correspond to an old break-up of the Nysa- 
Polana family. As a consequence, we can conclude that TC3 
comes from the Inner Main Belt, more particularly from the 
Background at low inclination or the Nysa-Polana family, 
and that these two sources could be genetically linked. 



6 POSSIBLE FORMATION SCENARIO OF 
2008 TC3 

Meteorite strewn fields with fragments of different miner- 
alogical types are very rare. Almahata Sitta and Kaidun 
llvanov et al.|[l983 ) are probably the only known specimen. 
In the scenario that mixing occurs by a coUisional pro- 
cess, the paucity of mineralogically-mixed meteorites sug- 
gests that the process that formed Almahata Sitta i.e. mix- 
ing the material between projectile and target is very rare; in 
most cases the proje ctile is pulve rized and leaves negligible 
traces in the target !|Meloshlll989l '). 

It is unclear which conditions allow for mixing between 
projectile and target. Impact velocity is probably a key pa- 
rameter. The average impact ve loci t y be tween asteroids in 
the main belt is 4.4-5.3 km/s (jBottke et al. 1994). If pro- 
jectile/target mixing were possible at these impact speeds, 
meteorites like Almahata Sitta would probably be frequent, 
which is not the case. Thus, we think that unusually low 
impact velocities are needed for mixing. This could prevent 



the target from pulverizing and could lead to macroscopic 
projectile fragments being implanted in the regolith of the 
impacted body or gravitationally bound to fragments of the 
target, if the latter is catastrophically disrupted by the im- 
pact. The fact that we give evidence that TC3 comes from 
the Nysa-Polana family (see Section Q or the Background 
at low inclination (see Section[5]), which are characterized by 
a mixing a taxonomic classes, also suggests that a specificity 
of the members of these families/regions may be unusually 
low collision velocities with projectiles that are also on low- 
inclination orbits. 

For all these reasons we did a systematic search for pro- 
jectiles that could hit Nysa-Polana family members or dis- 
persed asteroids of the Background at very low speeds. We 
did this search using the algorithm for the calculation of 
the intrinsic collision probability Pi between pairs of aster- 
oids, described in Wcthcrill (1967^. This algorithm, given 
the semimajor axis a, eccentricity e and inclination i of the 
two orbits (we use proper values for each pair of selected 
asteroids) assumes that the angles M, u, Q. (mean anomaly, 
argument of perihelion and longitude of the ascending node) 
of the two objects have a uniform probability distribution 
over the range — 27r; then it computes which fraction of 
these angles corresponds to the two objects being closer to 
each other than 1 km; finally, this fraction is translated into 
an intrinsic collision probability per year {Pi), using the or- 
bital periods of the two objects and assuming that they are 
not in resonance with each other. For our goals, we modi- 
fied this algorithm in order to take into account only orbital 
intersections corresponding to relative speeds smaller than 
0.5 km/s. Admittedly, this velocity threshold is arbitrary. 
Given the exceptional character of Almahata Sitta, we need 
a threshold much smaller than the typical impact velocity 
among random asteroids, that is to say a sub-sonic impact 
velocity, in order to preserve the target. 

Because most of Almahata Sitta fragments recovered 
are mainly made of ureilites, which we consider as analogs 
of B-class asteroi ds (as previously ment ioned in Section [1] 
and according to Ijennisken s et al.ll2010 !). we assume that 
the favorite scenario for the formation of TC3 involves a low 
velocity collision of asteroids of X and S classes (projectiles) 
with a B-class member (target) of the Nysa-Polana family. 

As for projectiles we considered all asteroids known 
in the ASTORB file (Asteroid Orbital Elements Database: 
ftp://ftp.lowell.edu/pub/elgb/astorb.html) of the whole 
Main Belt and we computed the probability of collision of 
each asteroid with a B-class object near the center of the 
B-class group of the Nysa-Polana family. Figure [6] shows 
the region of orbital element space where impacts can occur 
(impact probability >0). This figure shows that asteroids 
included in the region of low-velocity impacts have differ- 
ent SDSS colors and thus different taxonomies. As a con- 
sequence, collisions involving B-class, S-class, and X-class 
asteroids are likely to be possible. 

In the end, depending on the location of the target 
within the Polana group, we find that projectiles come from 
one of the following families: (i) Nysa-Polana, (ii) Flora or 
(iii) Massalia, located in the inner Main Belt, or (iv) Hestia 
which is in the central Main Belt, very close to the 3:1 mean 
motion resonance with Jupiter. More particularly, for targets 
with a semimajor axis (a) smaller or equal than 2.3 AU, the 
Flora family and objects of the background, very close to the 
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Figure 6. Distribution of proper orbital element for asteroids of 
the inner Main Belt: semimajor axes and eccentricities (top panel) 
or inclinations (bottom panel). For an asteroid target (black 
square), the region of orbital element space where impactors can 
come from is plotted in black. The location of the Nysa-Polana 
and Massalia families are respectively represented in blue and red. 
The color palette_ corresponds to the Parker's color distribution 
as mentioned in Parker ct al. and depends on the value 

of a* and i' - z' , with a* = 0.89(s' - r') + 0.45(r' - i') - 0.57, 
g' , r' , i' , z' corresponding to SDSS magnitudes at the following 
central wavelengths: 477 nm, 623 nm, 763 nm and 913 nm, respec- 
tively. Bluer dots (a* G [-0.25; 0] and i' - z' £ [-0.2; 0.2]) are 
related to dark objects while yellow to red colors (a* G [0; 0.25] 
and i' — z' £ [—0.2; 0.2]) gather brighter objects (see Figure|3]for 
more details). Data come from the SDSS Moving Objects Catalog 
4. 



Flora and Nysa-Polana family are the main source of pro- 
jectiles; for targets with 2.3 < a < 2.4 AU and eccentricities 
larger than ~0.15, projectiles in general are likely to come 
from the Nysa-Polana family; for targets with 2.3 < a < 
2.4 AU and eccentricities smaller than ~ 0.15, projectiles 
are likely to come from the Massalia family; for targets with 
a ^ 2.4 AU projectiles come from the Hestia family. Similar 
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Figure 7. Size frequency distribution (SFD) of the asteroid pop- 
ulation in the impact region. The SFD is plotted in terms of the 
cumulative number of asteroid as function of their absolute mag- 
nitude (H). 



results are found when studying collisions involving B-class 
asteroids of the Background at low inclination. 

Once identified the region of potential impactors, we 
proceded to estimate the likelihood and the frequency of 
impacts, which depend on the sizes, number density of as- 
teroids, and average impact probability in the region. The 
average impact probability in the region was derived by av- 
eraging the intrinsic impact probabilities (one of the output 
of our algorithm) of the asteroids in the region; we found 
Pi = 3. 10~^^ impacts/km'^/year. 

First of all, we estimated the number of impacts and 
the size of the impactors hitting TC3 during its coUisional 
lifetime assuming a diameter o f 4 m (given by its H value 
of 30.9 (| Jenniskens et al.l |2009| ) and assuming a geometric 
visible albedo of 0.05). The impact probability per projectile 
for TC3 is given by Pp — R^TcoiiPi with R the radius of TC3 
and Tcoii its coUisional lifetime estimated from iMarchi et al] 
(2006) {tcoU = 16.2 Myr). As a consequence, we find that 
TCa has an impact probability per projectile of the order of 
2. 10~^^ during its coUisional lifetime. 

To obtain the number of impacts on TC3 during its 
coUisional lifetime, we estimated the number of projectiles 
by calculating the size frequency distribution (SFD) of the 
population of the asteroids in the impact region and extrap- 
olating the corresponding function to asteroid sizes below 
the detection limit of current discovery surveys. Figure [7] 
shows the SFD in terms of the cumulative number of aster- 
oid as function of their absolute magnitude (H) . We used H 
values from the ASTORB file. In order to extrapolate the 
SFD asteroid sizes below the current detection limits, we 
calculated the equation of the straight line that best fits the 
base-10 logarithm of the known SFD and that is compati- 
ble with the size frequency d istribution of the whole Main 
Belt (see iBottke et al.|[20o3 ) (see Figure [T]). The line has 
equation y = aH + (3 where y is the \og-^Q{N{> H)) 10- 
base logarithm of the cumulative number of asteroids with 
H smaller or equal than a threshold value. The best-fit value 
for a and /9 are 0.5 and -4.5 respectively, implying about 10* 
asteroids with H < 17, 10^ objects with H < 21, 10^ with 
H < 27 and lO^'^ with H < 31, if the linear extrapolation is 
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assumed to be valid. We can deduce that for H < 31, we ob- 
tain only ~ 10~^ impacts with TC3. In the end, 1 impact at 
V ^ 0.5 km/s during TC3 coUisional lifetime would require 
1/Pp objects (i.e. 5. 10^'' objects). Using the SFD equation, 
this means we need a collision with an object of magnitude 
H ~ 42. This correspond to the size of a pebble of 2.4 cm of 
diameter. Clearly the contamination brought by low speed 
impacts is very small (the volume ratio between TC3 and 
the pebble is 1.1 10~^) and certainly not at the level of 20- 
30% as implied by the number of chondrites found amongst 
the ureilites of the Almahata Sitta meteorites. 

Much more likely, TC3 was a fragment of a larger aster- 
oid that was liberated by an impact near the z/g resonance 
few millions years ago. In fact, we note that (1) the dy- 
namical time to deliver asteroids from the resonance to 
orbits similar to tha t of TC3 is of the order of 19 Myr (see 
ICampins eraLll2010l . for 1999 RQie coming from the Polana 
region); (2) the cosmic ray ex posure age of Alma hata Sitta 
meteorites is of 19.5±2.5 Myr (jWelten et al.ll20l"ol 'l: as a con- 
sequence, TC3 spent most of its life as a single asteroid trav- 
eling from the ua resonance to the orbit that brought it to 
impact the Earth. 

It makes sense to investigate what was the contamina- 
tion brought by low velocity impacts on the parent body of 
TC3. The average impact probability implies that we have 
~ 1.4 10~^ impacts/km^ over the age of the Solar System, 
which correspond to 1 impact of an asteroid of ~ 27 tak- 
ing into account the SFD of the objects in the region. The 
diameter of an asteroid with H=27 is between 10 and 25 m 
depending on the albedo. Again the level of contamination 
is very little, of the order of 1. 10~^. 

When we take into account the necessity of mixing me- 
teorites of at least three different mineralogies, we conclude 
that the formation of TC3 by means of low- velocity collisions 
is not a realistic scenario. 



7 CONCLUSION 

From our study, it seems that the Nysa-Polana family and 
the Background at low inclination are good candidates for 
the origin of 2008 TC3 and Almahata Sitta. First of all, 
as mentioned in Section [51 the Nysa-Polana family is lo- 
cated close to the secular resonance which is the favorite 
route leading to primitive NEOs and more particularly to 
asteroid 2008 TC3. Moreover, the proper inclination of the 
Nysa-Polana family is very similar to that of TC3, which 
should have been maintained during the transfer of TC3 
to the NEOs region through the i/a secular resonance. We 
also know, from our algorithm of spectral classification (Sec- 
tion [S]), that the Nysa-Polana family gathers the 3 spectral 
classes (S, B, and X), which are proposed analogs to Alma- 
hata Sitta fragments (Section [4]) under the hypothesis that 
ureilites are linked to B-class asteroids. More specifically, (1) 
the mean spectrum of B-class asteroids of the Nysa-Polana 
family is spectrally matched - at least in the visible - with 
available spectra of ureilitic fragments of Almahata Sitta, 
(2) considering space- weathering effects, the mean spectrum 
of S-class asteroids of the Nysa-Polana family is compatible 
with the spectrum of the H-chondrite fragment #25, (3) 
a good agreement is found between X-class asteroids and 
enstatite chondrites from other meteorite falls (we remind 



that enstatite chondrites are part of Almahata Sitta). Of 
course, a comparison with enstatite chondrites from Alma- 
hata Sitta fragments would be very useful to get a definitive 
match between the Nysa-Polana family and Almahata Sitta. 
In Section [5] the same kind of work was performed for ob- 
jects of the inner Main Belt coming from the Background at 
low inclination (i < 8°). We concluded that these dispersed 
asteroids could also be at the origin of TC3. 

In Section [6] we tried to explain the formation of TC3 
by low velocity impacts (below 0.5 km/s) between different 
mineralogies in the neighborhood of and within the Nysa- 
Polana family. Selecting TC3 as the target asteroid (d=4 
m), we find a probability per projectile about 10~^^ im- 
pacts during its coUisional lifetime (i.e. in 16.2 Myr). As a 
consequence, impacts at low velocity are extremely rare and 
there is little chance that TC3 was formed by low-velocity 
impacts in the current asteroid belt. This implies that the 
heterogeneous composition of the parent body of TC3 has 
to be inherited from a time when the asteroid belt was in a 
different dynamical state, most likely in the very early Solar 
System. One could think that an asteroid of ureilite com- 
position was contaminated by impactors of different nature 
when the asteroid belt was still massive and dynamically 
cold, so that mutual collisions were frequent and at low ve- 
locity. However, this view is probably simplistic. In fact, a 
body of ureilite composition needs to be formed in the inte- 
rior of a large carbonaceou s asteroid which underwe nt sig- 
nificant thermal alteration (|Singletarv fc Grove|[20o3 ). This 
asteroid needs to have undergone a coUisional disruption to 
expose the ureilite material in space. The same i s true for the 
bodies of Hn composition, with n larger than 3 (|Gopel et al.l 
[1993). But coUisional disruptions require large relative ve- 
locities, in contrast with the view of a dynamically cold belt. 
Thus, the asteroid belt could not be overall dynamically cold 
when the parent body of TC3 formed. 

These considerations suggest that, conversely to what 
is usually thought, accretion and coUisional erosion had to 
co-exist for some time in the asteroid belt. For this to be pos- 
sible, presumably there had to be still a significant amount 
of gas in the system so that, although large asteroids could 
be on dynamically excited orbits, the orbit of their small 
fragments were rapidly re-circularized by gas-drag. Conse- 
quently, the mutual relative velocities of these fragments 
were small and a new phase of accretion was possible for 
them. 

We remark that the heterogeneity of TC3 is not at the 
microscopic level; each of the meteorites delivered to the 
ground are of a distinct class. Thus, TC3 seems to be an 
agglomeration of meteorite-sized (i.e. few dm) pebbles of 
different nature. Pebbles of this size are strongly coupled 
with the gas and are extremely sensitive to pressure gra- 
dients. They play the key role in the new models of plan- 
etesimal formations, based on the concentration of dm-size 
pebbles in the eddies of a turbulent disk and on the process 
of streaming instability (|johansen et al.ll2007l . [200^ ). These 
models of planetesimal formation in a turbulent disk seem 
a priori to be particularly favorable to explain the coexis- 
tence of coUisional erosion and accretion. Large planetesi- 
mals are formed by the concentration of a large number of 
pebbles, forming self-gravitating clumps. Once formed, the 
orbits of these large planetesimals are rapidly excited by the 
stochastic gravitational perturbations exerted by the tur- 
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bulent disk (|lda et alJliooi : iMorbidelli et al.l 120091 '). If the 
threshold for colHsional break-up is achieved, the pebble- 
size fragments of these large bodies are re-injected into 
the game: by being concentrated into new eddies they can 
give origin to new large planetesimals and so forth. Admit- 
tedly, quantitative work is needed to support this scenario; 
also, other more classical plan etesimal formation mecha- 
nism in presence of gas drag ( Wcthcrill fc S teward \W9^ 
iKenvon fc BromlevI 12004 IWeidensc hilling 20li|) might ex- 
plain the coexistence of erosion and accretion as well. 

In this respect, it will be important to understand from 
the observational point of view if macroscopic heterogene- 
ity as that of TC3 is the exception or the rule among as- 
teroids. TC3 is the first object of this kind that has ever 
been observed, but it is also the first fall of an asteroid on 
Earth documented live and for which an extensive and ex- 
haustive search for meteorites has been done. So, it might 
not be as rare as one could be tempted to b elieve. Indeed, 
a sec ond similar case has just been reported (|Spurnv et al.l 
|2012| ). Now that the possibility for macroscopic heterogene- 
ity is recognized, careful investigations (also conducted by 
remote sensing techniques) may reveal additional interesting 
cases. Understanding which fraction of the asteroids are of 
primary or secondary accretion will be a fundamental step 
to constrain the asteroid formation and evolution models. 
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